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I. INTRODUCTION 



As one derives the effective locality property of the fermionic QCD amplitudes [U [2], 
a bizarre factor of 5^(b) comes about, in an unexpected place, where b is the impact 
parameter of two scattering quarks. Inspection shows that this curious factor is not related 
to any systematic approximation scheme, such as eikonal or quenching. Rather, it shows up 
as an unavoidable feature of any treatment, exact or approximate, and is likely to suggest 
the neglect of transverse fluctuations of quarks confined inside hadrons. In this respect, 
quite interesting perspectives offer themselves [2]. 

In this article, our aim is twofold. In the first place, it is made clear that this 5^ 2 \b) can 
not be avoided, and a rigorous mathematical proof is proposed to support this claim. Then, 
as a second step, in eikonal and quenched approximations, it is shown how the powerful 
Random Matrix calculus can be used to obtain the generic, closed form expression of QCD 
fermionic amplitudes in the strong coupling regime. Eventually, the forms so obtained will 
be used to revisit and improve upon a previous estimate of the quark/anti-quark binding 
potential and related aspects [3]. 

Quenching and the eikonal approximation will be used throughout. Both seem quite 
appropriate in the context of an energetic scattering process: The former because one may 
think that the vacuum back-reaction has not enough time to influence the process; and the 
latter, which can be defined in a Lorentz invariant way, because it preserves unitarity, a 
most important requirement of a non-perturbative analysis. 



II. THE CONSTRAINT OF DELTA-FUNCTION 

Using standard functional notations and manipulations pQ, and small re-arrangement 
that guarantees rigorous gauge invariance for all scattering quark amplitudes [2], the QCD 
generating functional can be written as 

Z\j,r},fj] = M ei f&P* J d[ X ] e^* 2 e SA e^ /xf+a /ad^a 

(1) 



x e 



JfjGc[A]r, L[A] 



where M is a normalization constant, and the and r] a , fjp stand for bosonic and fermionic 
sources, respectively; and the causal gluonic propagator D c may differ by a factor of 2i from 
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the more customary Feynman prescription. The 'Linkage Operator' 

D A = -'- [ d 4 xd 4 y— ^— rtP\ ab (x-y)— ^— (2) 
involves the covariant gluonic propagator, with gauge parameter £ 



BP\ ab =5 ab (-d 2 )- 1 



c 



g»u - (i - C) 



Q2 



(3) 



The X- held appearing in Z[j, r], rj\ is the Halpern variable [5] used to linearize the exponential 
of the original F 2 -dependence that belongs to the original QCD Lagrangian density: 

e-i / p2 =Af x J d\xU d S^+i /** "*Jk (4) 

Proceeding in this way, the A^-gauge field dependence are made gaussian instead of cubic 
and quartic, and the linkage operations can be carried out exactly 

In Eq. 0, the original fermionic fields have been integrated out: This has given rise to the 
term involving the fermionic propagator G C [A] in the gauge field background configuration 
A, as well as to the logarithm of the fermionic determinant, LL4], that is, 

G c (x, y \A) = (x\[m + Y % - % g A; T a )Y l \y) , (5) 

and, 

L[A] = Trln[l -ig( 7 AT) S c ], S c = G c [0]. (6) 

Note that the property of effective locality is not discernable from the generating functional 
itself, but only after the sums over all gluon exchanges have been performed between those 
quark lines that define the amplitude under consideration; that is, on the full set of 2n-point 
fermionic Green's functions. Though equivalent, and thus not a crucial difference, a form 
of the effective action is under investigation, on which effective locality could be read off 
directly. This point would be important when this formulation of QCD amplitudes generates 
a possible most dual QCD formulation [I], or in relation to the most recent results of Ref. [?]. 

In order to be able to display the property of effective locality, one must devise convenient 
enough representations for the functionals G C [A] and L[A]. Schwinger-Fradkin representa- 
tions can be used; for example, a ' mixed' (configuration and momentum space) expression 
of G C [A] is given by 
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/■oo 

(p\G c [A]\y) = e~ ip - y % \ dse~ ism2 e ^ Trln(2/l) 

Jo 

x J d[u]{m - z 7 • [p - gA(y - u(s))}} e*£ ds ' [u ' (s ' )]2 e ip ' n(s) 

x ^ e g/ s d S 'a-F(2/- U ( S ')) e -i 3 / s d S '«'( S ')-A(y- U ( S '))^ ? ^ 

and similarly for LL4] [2]. In the representation above, one has 

d d 

h(st, s 2 ) = SiO(s 2 - si) + s 2 Q(s 1 - s 2 ) , h (s 1 ,s 2 ) = ——6(s 1 -s 2 ), (8) 

OSi os 2 

where a^ v = |[7 M ,7^]- This expression ends up with a subscript '+' to mean that time 
ordering with respect to s' is in order, because of the S'L r (iV r c )-Lie-algebra valuation of the 
gauge potentials A^. 

The Fradkin field variables satisfy 

w M (s)= / ds'v^s'), u„(0) = 0. (9) 
Jo 

On the mass shell, v M (s') can be interpreted as the particle 4-velocity p fl (s)/m. The relevant 
eikonal-approximated form of (J7]) that will be used here is 

(p\Gc[A]\y) 



x 



/■oo 

/ ds e~^ m2+ ^ {m - iY [p, - 9 A a (y - sp)\ a ] } 
Jo 

exp |s jT ds'a^ F« v (y - s'p)A"} • exp {+2^ jT ds' A^y - ^)A a ] J J (10) 



The full derivation of the effective locality property will not be repeated here since its 
complete non-approximate proof can be found in Ref. [2]. Rather, focus will be put on 
specific and important technical aspects. 

In a first attempt to the EL property of QCD [lj, one finds an exponential factor of 



% 

2g 



d'xQl{x){f. X {x)YJ2 QUx) (11) 



representing, since all gluon interactions have been summed upon, the action of an effective 
local interaction term between two colored currents. The overall currents are 

Ql{x) = -d v xlM + g[RU x ) + ^O*)]. (12) 
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where the Rf , i = 1,2, are the associated currents of quarks. Without loss of generality for 
our purpose, in the following, quark species will be reduced to one flavor only. 



In a non-perturbative strong coupling regime, g 3> 1, the leading part of (11) reads 
accordingly [TJ [2] 



l -g J d 4 w J dsi J ds 2 u'^{s 1 )u' u {s 2 ) (13) 

x 5^(w - Vl + u( Sl )) S^(w -y 2 + u(s 2 )) , 
where the O a (si) and f2 fe (s 2 )-variables are used so as to extract the A^-field dependence out 



of the ordered exponential of (10) [2]. This expression is obtained in the approximation of 
quenching. However, it is important to emphasize that the full non-approximate expres- 
sion [2j would manifest exactly the same technical intricacy as the one under consideration, 
and this is why the point can be made using this simplified example. 



The technically important aspect of (10) is the following: What is to be done with 



5 (4) K - vi + ^\yi -V2 + u(s 2 ) - u( Sl ))i (14) 

That is, basically, how should one interpret such a factor as 5^\C st + u(s 2 ) — u(si)) ? 

At face value, for any given pair of values (si, s 2 ) G ]0, s] x]0, s], and any pair of arbitrary 
functions (u, u), each belonging to some infinite dimensional functional space, the probability 
of coincidence of u(s\) with u(s 2 ) + C st is likely to be infinitesimally small, if not identically 
zero. 

Besides, in the 2-body center of mass system (CMS), a somewhat heuristic manipulation 
of the 0- and 3- delta constraints would suggest that C st = 0, and write [2] 

6(s }S(s ) 

5(u (s 2 ) - u (si)) S(u 3 (s 2 ) - u 3 (si)) = + • • • , (15) 

where the dots are to be evaluated shortly. It is assumed that the u, u Fradkin fields are 



C 1 (]0, s], ]0, s] — > 1R 4 ). Out of 5^(u(s 2 )— u(si)), and in view of (15), this leads to a remaining 
constraint of 

S {2 \y 1± -y 2± )^5^(b), (16) 

— # 

where b is the impact parameter, or transverse distance between two scattering quarks. Given 
the very place where it comes about, as a multiplicative factor in an exponential, this S^(b) 
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is awkward. In Ref. [JJ, it has been suggested that this factor is a remnant of an implicit 
existence of asymptotic quark states, an assumption which, beyond the stage of perturbation 
theory, cannot be maintained in QCD [5], neither theoretically nor experimentally. 

This point of view leads to a change of 8^ 2 \b) into a modified impact parameter distri- 
bution of form [2], 

8®(b) -+ v (b) = £ e -0*)« E 9 e « 1, (17) 

where \l is a typical mass term used to calibrate the transverse momenta distribution of 
quarks inside a given bound state. In the transverse plane, bound quarks are endowed with 
transverse momenta taken as independent random variables, exponentially suppressed above 
a given mass parameter /i. That is, describing the non-perturbative regime of QCD, in order 
to avoid the nonsensical 8^ 2 \b), it is necessary to introduce a mass scale where there wasn't 
any before, in agrement with the considerations of Ref. First fits, with respect to a 
model pion QQ and a model nucleon QQQ [2], indicate a value of \i close to the pion mass, 
and a value of £, y/2/16, on the order of 0.1. 

Eventually, starting from quark propagation as ordinarily conceived, a la G c (x,y\A), 
in the non-perturbative bound context, one could be lead to think of transverse quark 
propagation in terms of Levy-flight, a picture which, in view of color confinement, may look 
rather appealing indeed (Levy flight distributions are stable probability distributions. They 
include and generalize gaussian distributions for independent random variables, and comply 
with a generalized 'central limit theorem' of statistical physics [TO]). 



Whatever its interpretation, though, the latter will be ruined if ( 15 ) does not provide a 



reliable enough evaluation of (14). In order to deal with that issue, one may skip to the most 
achieved realization of a functional space, that is to the Wiener functional space [UJ. Note 
that even though that space may be not the only one to be conceived (and/or taken into 
account, according to Ref. [7]), there is no restriction at all in working out a proof within 
that space, because Fradkin vectorial fields are real valued. Then, the following theorem 
can be proven: 

Theorem: For all pairs (si,S2) £ ]0,s]x]0, s], 

m®m ({(u,u) G C 0,s x C °' s " | u( Sl ) = u(s 2 )}) = 0, (18) 
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({(u,u)eC° ' s xC? 







u(0) = u(0) = 0}) = 1 



(19) 



with m, the Wiener measure on Cq' s , whereas m <g> m is taken as the Wiener measure on the 
product of spaces Cq' s x Cq' s , endowed with the topology product. This is made possible 
thanks to the independence of the random variables u and u that are here taken to represent 
either of the 2 possibilities uo,3 and ^0,3. 

This theorem proves that the right hand side of ( [i~5"| is necessarily proportional to 
S(s\)5(s 2 ). Dimensional and symmetry arguments (the right hand side must be symmetric 
both in the exchange of indices f> 3, and under the combined exchange Si s 2 and 
u -H- u, whereas the normalization is fixed by comparison to the eikonal approximation pQ) 



can be used to complete the right hand side of (|15|). One finds 

S(u (s 2 ) - U Q (S 1 ))S(U 3 (S2) - « 3 (si)) = I 



S(s 1 )S(s 2 ) 



u' 3 (si)| \u' (s 2 ) 



+ 



u' ( Sl )\ \u' 3 (s 2 ) 



(20) 



and because of the CMS-relation u' (si)u' 3 (s 2 ) = u' 3 (si)u' (s 2 ) , which holds true at eikonal 



approximation, (20) gets reduced to (15), without dots. 
The proof is as follows: 

Let A be the set {(u,u) G C°' s x C°' s 



UlSi, 



m(s 2 )}. One has A = H^Li Ai> where 
A n = |(«,«) e C U ' S x C U ' 5 I - ^ < -«(s 2 ) < +ij (21) 
Because of the obvious inclusion, Vn, A n+ i C A n , one can write 

m <S> m(A) = lim m <g> m(A n ). 

n— >oo 

Now, X n = m®m(A n ) is given by 



(22) 



+00 



00 
+00 



| (w, u 


) c c °' s x cr 


1 u ( s l) 


dx 


rx+- 

^ n 


dy 


a/27TSi 


a/27TS 2 


dx 


e 2 ^ f n {x). 




\f2-KSi 





- < U( S 2> < U{S1) + - 

n 



n 

,-y 2 /2s 2 



(23) 



(24) 



Since one has the obvious inequality of 

1 



e 



< 



e 



(25) 
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one gets immediately 

^ = (26) 

as a mere consequence of the Dominated convergence theorem, allowing one to take the 
n = oo-limit under the integral, whereas the second equation follows from normalization. 

The presence of this surprising factor of (b) is thus unavoidable, and that issue bounces 
back and forth with the related emergence of the £ parameter: Again, if this parameter is 
really an unavoidable output of this approach, more work is in order so as to make explicit 
the relations of £ appearing in the confining potential between a pair of quarks separated 
by a distance r, to wit |2J 

V(r) ~ (27) 

to the Levy-flight propagation modes of confined quarks, as well as to an unexpected non- 
commutative geometrical aspect of the scattering transverse planes. 



III. FERMIONIC QCD AMPLITUDES IN EIKONAL AND QUENCHING AP- 
PROXIMATIONS 



As a result of EL, one finds a 2-body scattering amplitude [TJ |2] proportional to 



d n ai / d n a 2 I d™^ / d n n 2 e- iai - ni e~ ia2 - n2 e iai - T e M 



Ka<n 



N II / d tC»] det[gf- X (w)] * < 

0<H<i/<3 



(28) 



where n is the shorthand for n = N^—l, and the • is used to mean an Euclidean scalar product 
in a relevant adjoint representation space: For example, the product [/ • x] bc — Y2 a / 
The T a are the n hermitian traceless generators of the S , 6 r (A r c )-Lie-algebra, taken in the 
n-dimensional adjoint representation of SU(N C ). In this representation the T's are given by 
the structure constants of the S'6 r (A^ c )-Lie-algebra: (T a ) bc = —ifabc = —if abc . 



In the second line of (28), the last exponential argument is nothing else than (13) evalu- 



ated with the help of identity (|20|). In eikonal approximation, one gets in effect 

1 



exp 



W6)<(0)<(0) 



lit 



[/•x(«0r|£n;(o)n$(o) 



(29) 



where u'(0) = P\jm and u'(0) = P-ifm, and where, as advertised, the original 5™ (ft) has 



been turned into the distribution (pip) of (17). In the CMS of the 2 scattering quarks of 



momenta Pi and P2, (29) reduces to the last argument in the exponential term ending the 



expression (28) 



It is worth recalling that the peculiar form of (28) is non-trivial and originates in the 



property of Effective Locality: Having selected w = (0,y±,0)cMS with y± = y\± = —ij2± = 
^b, as the only spacetime point where the interaction takes place, the contributions attached 
to all of the other points are eliminated. The result is that there is no spacetime integration 



in the argument of the exponential term at the end of (28). And accordingly, there is no 
spacetime integration either, for the first term. 

This can be understood in a straight forward way by recalling the standard procedure of 
a generating functional construction [5] , where the spacetime manifold is broken up into an 
infinite series of infinitesimal cells. Since the interaction, of the contact type, takes place at 
a unique selected point, w, the contributions of all of the other infinitesimal cells, centered 
at different points, just contribute an infinite product of multiplicative factors of 1, thanks 
to the normalization of each cell contribution. 



However, the final result for the amplitude will depend on 5 (cf. (31)), and hence a 
non-trivial result will depend upon forming a connection between the vanishingly-small, 
mathematical 5, and a physically meaningful 5 p h consistent with experimental measurement. 
That this is both necessary and can be achieved follows from Effective Locality (EL), and 
can be seen as follows. 

After EL has forced only a single xe to appear in the interaction exponential, with all the 
other xi'i V 7^ ^ integrals trivially evaluated and combined with their normalization factors 
to give a net contribution of 1, we are left with the integral 

Af(8) J d n X i e^ x ' (30) 

over the interaction part of the integral which, when evaluated, depends upon 5. Now define 
a re-scaled variable xt by the relation 

S 4 ph X 2 e = S 4 xl (31) 
so that the effective normalization and the new measure combine to give 



m) [ -f] I d n Xi e^«. (32) 



But AT (6) = C st ■ 5 2n , so that Af(5) (5 ph /5y n = Af(5 ph ) as expected. 



EL forces us to think in physical terms, and the only remaining question is how to choose 
<5 p h, which quantity will now appear in the interaction integrand (multiplying <p(b)); and 
we let Quantum Mechanics provide the answer: That contribution corresponding to a time 
separation should be chosen as 1/E, that corresponding to a high-energy CM longitudinal 
coordinate should be 1/pl — while that corresponding to each of the transverse coor- 
dinates should be where \i is on the order of the pion mass; and hence S^ h = (1/E/i) 2 . 

To paraphrase Schwinger [6], just as reorganization is needed, in principle, to make the 
connection between interacting fields and asymptotic Abelian fields, so an effective renormal- 
ization is required to pass from non-Abelian fields, whose quanta do not have asymptotic 
states, to a practical description in which their scattering or bound states do have a full 
asymptotic description. That effective renormalization is defined by the replacement of the 
mathematical 5 by the physical 5 p h- We subsequently suppress the subscript of <5 p h, and the 
bar notation of % 



As it stands, though, the expression (29) confronts us with a puzzling issue, namely, the 
inverse of the quantity [/ • xi w )]- Taken as some relevant matrix, in effect, it is immediate 
to realize that zero belongs to the spectrum of [/ • xi w )}, so that there seems to be no way 
out of this dead end. However, two most favorable circumstances allow one not only to cope 
with that issue, but also to reach the level of actual sensible calculations. 

The first one of these two circumstances is as follows. In order to: 

1. Invert, at least formally (for the time being) the product [/ • x( w )]> 

2. and take into account the contributions of all of the 6n- -fields entering [/ • x( w )] 
at the spacetime dimension of D = 4, 

3. as well as the full symmetry of the product [/ • x( w )] under the combined exchanges 
a <->■ b, // <H- v [I], 

it is appropriate to think of (f -x) as a tensorial product Yla=i X a>1 v®T a of operators acting 
on the product space {//, u} £g> {a, b}. The first space, {/x, v}, is the vectorial Minkowski 
space (as opposed to the affine Minkowski space, the space of event-points, to which the 
former space is tangent), and the second, {a, b}, a S77(iV c )-adjomt representation space with 
dimension of n = N% — 1. Their product is a Dn = iV-dimensional vectorial space that can 
be endowed with the product basis {(fx, (p 2 , ■ ■ ■ , <^jv} constructed out of the bases of {/i, u} 
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and {a, b} spaces, on which \ ail v and T a act as linear operators. Then, in this product 
space, and at real-valued Halpern-fields (a generalization to complex-valued fields has been 
proposed recently [7]), the operator X]"=i X afJ- u ® T a is represented by a A^ x A" matrix equal 
to a real anti-symmetric matrix M multiplied by i, that is, a traceless hermitian matrix, 

n 

X a \ ®T a = iM, M i3 = -Mji (33) 

a=l 

As will be noted shortly, this result has an important consequence for calculation purposes. 
In Minkowski spacetime, which is pseudo-Euclidean, it depends in a crucial way on the 
fact that the original antisymmetric Halpern fields, are taken as {/i, z/}-space endomor- 
phisms, that is under the very form of mixed tensors, x a>1 v - 

The operator x®T acts upon the A^-component vectors w'(0)(g>f2i(0) = V\, w'(0)<g)f22(0) = 
V2 of the product space, and in the basis {ipi, <f2, • • • , <Pn}, these vectors have components 

v x = (u' nl,u' nl . . . , u'^f, v 2 = (u' nl,v/ nl . . . , u' 3 n^) T , (34) 

where, for short, the dependence on Si = 0, i = 1,2, are hereafter omitted. Note that like 
the fij's, the V^s are also dimensionless objects since, as quoted before, u'(0) — P\/m and 
u'(0) = P 2 /m. 

There therefore exists a unitary matrix U that takes iM to its diagonal form 

Wt(zM)W = diag (...,&,...), i = 1,2,..., N, & e Sp(zM), (35) 
where is the hermitian conjugate of U, and likewise M _1 to its diagonal form also 

W t (zM)- 1 W = diag(...,^,...), z = l,2,...,AT. (36) 

Si 

It is easy to check that the spectrum of (iM) contains the eigenvalue £j = at multiplicity 



greater or equal to 1 so that, for the time being, (36) can only be given a formal meaning. 



This allows one to re-write ( 29 ) as 



exp (-Mb)^ VI \{iM)-% vj} , (37) 
and then, by redefining the V^-fields, 

Vi^Vf = UVi, (38) 
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one gets for (37) the expression, 



m 



N 



exp | -9<p(b)-^ 



At eikonal approximation, fixed u', u' = Pi/m, Pijm, one has d n fli 



(39) 



m 



-n a n 



d N (Pi®Qi), and 



likewise, the Jacobian of the transformations (38) is trivial. By using the reformulation (33), 
the integration with measure 



l<a<n 

n 

0</i<i/<3 



(40) 



which appears in (28), gets translated into an integration over the algebra of hermitian NxN 



matrices, with measure [13 

n 

d(E X - * - ® T °) 



a=l 



d(iAf ) = d(iM)nd(zM)i2 • • • d(z'M) jviv 
<9((iM) n ,--- ,(iM)iw) 



<9(£i,--- )Cjv,Pi,-- - ,Pn(n-i) 

N 

n d & n & - & r - d^(v-i) /(p) 



d£i • • • d^AT dpi • • • dpAT(Ar_i) 



(41) 



1=1 



l<i<j<N 

The iV eigenvalues of (iM) are real, and f(p) is an (unspecified) function of the N(N — 1) 
extra parameters which complete the parametrization of a given hermitian matrix (iM). 
Mathematically, this passage from an infinite dimensional functional space of integration, to 
a finite dimensional matrix space is non-trivial: It is a direct consequence of the Effective 
Locality and can be proven by relying on the measure image theorem [TT) IT2"] . Since nothing 
in the integrand depends upon the parameters pj, their integrated contribution factors out 
and disappears with the normalization so as to leave an integration over Sp(iM) G M only. 



The important consequence alluded to after (33) shows up here, in (41), with that power 
of k being equal to 2 in the case of an hermitian matrix iM, changing the unpleasant absolute 
values of |£j — into the analytic functions — £j) 2 . 



The second line of (|28|) can accordingly be re-written as 

o 8iV r S ' 
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N 



n ^ 



i=i 



n 

X<i<j<N 



e -«K«^pE ( ^ f (42) 



where, for the exponential term, the relation XnuX allL ' = 2iV c Tr (iM) 2 has been used, and 



the same symbol as in (28) has been retained for the normalization constant. Keeping in 
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line with [13], the normalization constant may be so defined as to normalize (42) to unity 
when the whole set of functions 



1 

e ^ 



(43) 



is replaced by the set of constant unit functions, /, = 1. In order to allow for a concise 
and closed form expression, one can take advantage of the method used in [13]; a simple 



re-scaling of the eigenvalues is necessary, that is, £| = yl/8N c £i, taking (42) into the form 

N 



_N_ 

9 2 



(44) 



in which one has a Vandermonde determinant of 



n (e t -^) 2 = Poi(6,--.,^) 

l<i<j<N 



(45) 



with Pol(£ 1; . . . ,£n) a sum of 3^^ 1 - )//2 monomials in the £j's, of overall degree N(N — 1) 
each. Introducing the normalized harmonic oscillator wave functions, 

, , 1 



e-^ x2 H k (x) 



(46) 



7T 



where Hk{x) is the Hermite polynomial of order k, the expression 

^ 2 (6,..-,^):=^ 2 e-^^ 2 J] te-O) 2 , 

l<i<j<N 

can be written as 

N-l 

Pn 2 (£l> • • • j 6v) = CjV 2 ] J 



(47) 



fe=0 



2 fc A;! 



det(^_!(Vi^))] , j,k = l,...,N,(4S) 



where the constant C^ 2 is defined so that the integration of (47) over all of the N eigenvalues 



£i yields unity. Since P/v 2 (6i> • • • >6v) is nothing else than the integrand of (42) modulo the 



prescription (43), this would identify the normalization constant of (44), A/", as just this 



constant Cn 2 , which reads [13J 

C-i = 2 -( 1 /W-i) (iV)!7r f Jjjl 



N-l 



(49) 



3=0 



However, contrary to the standard expression of [13J, (47) displays a factor of i in its 



exponential term, which, in (48), gets reflected in the \fi ^-dependence of the harmonic 
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oscillator wave functions. Relying on Cauchy's Theorem, though, it is elementary to prove 



that the integration formulae of (47) and (48) defined without y/i, can be analytically con- 
tinued from real valued £j to the complex range of z = £,je l6 , < 9 < tt/4 (see the Appendix 
where a similar procedure is followed). In particular, using this analytical continuation it 
is possible to check that the harmonic oscillator wave function normalization undergoes a 
simple re-scaling of 



+oo 



dx (p k (x)(pj(x) = 5 k j 



which, proceeding along the standard steps of 



dx (pk{xVi)(pj(x\A)) 



1 



>hji 



(50) 



allows one to conclude that the normal- 



iV-l 



ization constant of (44) is given by 
and therefore, 

Pn 2 (€i,- ■■ ,6v) = T^y {det ^fc-i(Vi^))} , j,k = l, 



(51) 



N. 



(52) 



Eventually, by inserting (52) into (44) it is easy to obtain the 4-point fermionic amplitude 



in the most elegant and concise closed form of 

ff -?(8iV c )*-«^{det[F^>)]} 2 , 
where 'the matrix elements' are given by 



z = 0,l, 



Fi(B 



0> 



• N-l, j,= l, 



0> 



and the function v, by 



V8N c Ep 



(53) 



(54) 



(55) 



It is remarkable that exactly the same equations as Eqs. (53), (54) and (55) express the 



closed form of any fermionic amplitude: In the general case of a fermionic 2m-point function 
in effect, all of the previous steps can be followed through, and the only change turns out to 
be the redefinition of the constants B^K Within obvious notations (\pk\ = \pi\ '■= Pkh an d 
E k = Ei := Eki), 

,2 



B 



(./) 

2 m 



9 



m 



Kk<l<m 



Em Ph 



v'lv'l 



(56) 
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where again, the restriction k < I in the summation, discards the self energy effects that 
will be dealt with in another article. In Ref. [20], though, a glance at this issue is obtained: 
Amazingly, the same ultraviolet renormalization structure as in the usual QCD perturbative 
regime comes about, and it is interesting to note that this point is in line with the explanation 
proposed in Ref. [18], based on the AdS/QCD correspondence. 

Completing the first most favourable circumstance that was encountered at the level of 



Eq. (33), the second one is met in Section IV, where, illustrated in the simplest case of 
m = 2, it is shown that the matrix elements Fi(B^>) of Eq. (54) exist, are calculable, and 
their further auxiliary dependence integrable also. 

For all i, j indices, in effect, the matrix elements Fi(B^) turn out to be given by linear 
combinations of Meijer's special functions, elsewhere referred to as Fox's //-functions [TJ. 
To wit, by anticipating a bit of the next Section, 



Ft(B 



0> 



1 



JS(t/2) 

2~7^ 



■l) k (2k- l)!!Cf y/-8iN t 



(57) 



^03 



i(B 



4A; + 3 1 



.0 



where E(i/2) is the integer part of i/2, and the Cf are the binomial coefficients 



At a formal/algorithmic level, Eqs. (53)-(57) make up an important result of this article. 



Extended to the more complex case of eikonal and quenched QCD, Eqs. (53)-(57) seem to 
comply with a general conjecture which, in Ref. [7], is illustrated in the cases of simpler 
scalar models only. 



IV. G™ n -MEIJER'S SPECIAL FUNCTIONS 



"pq 



The general and closed form expressions given above may be most appropriate to numer- 
ical simulations. In this Section though, for the sake of illustration and a further speculation 



on non-perturbative QCD, it is enough to get back to the expression (42), here reproduced 
up to its overall multiplicative constants: 



+oo 



N 



n » 

X<i<j<N 



(58) 



The Vandermonde determinant being a polynomial in the eigenvalues, one has to cope with 
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generic integrations of the type 



(59) 



where ft G IN. Defining 



da; x p e x2 * 



(60) 



for p > and b > 0, one has that (60) evaluates to [I 

'6 2 



Jp(&) 



1 



4 



p + 1 1 



.0 



(61) 



where G^(x\bi, 62, ^3) is a Meijer's special function, and the second of the two most favorable 
circumstances alluded to above is the following: 



Whereas (61) holds true at values of p and b such as specified above, it turns out that 



Gl®(x\bi, 62, ^3) is analytic in x, and that the parameters appearing after x can be continued 
to complex values [Ej. 

On the basis of these analyticity properties, it is therefore possible to give the generic 



integrations of (59) the meaning of 



d£i „ » *2 „,„ V 'l V '2 



(62) 



/o30 
^03 



i(B 



(i)\2 



2ft + 1 1 



4 '2' '' 



2^ UJ V 16 

where the equality, demonstrated in the Appendix, is a mere consequence of Cauchy 's The- 



orem. 



and the constant is defined in (55) 



Now, each of the is on the order of g<p(b), and the small g<p(b) limit can therefore be 
investigated on the basis of the small z expansion, \z\ <C 1, 



4 u'Unb 3 -b h ) 



z bh x F 2 (b 1 ,b 2 ,b 3 ;-z), 



h=l 



+ b h 



(63) 



where the prime on top of the product indicates that the pole value bj = b h is omitted, and 
likewise the asterisk in the 0-^2 generalized hypergeometric series, indicates also that the 
value bj = bh is omitted: 



)F 2 (1 + b h - 61, ...,*, 



l + b h - 6 3 ; -z). 



(64) 
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For all bh, one has, at \z\ <C 1, 



)F 2 (1 + 6 h - 6 ls . 



l + b h -b 3 ;z) = l + 0{z). 



(65) 



For the situation considered, (61), the parameter sets (61,62,^3) always read as 



( 2g 4 +1 ; \ i 0); so that i n the small z-limit, the Meijer's special function (61) yields, for ex- 
ample, 

'(fl w ) 2 .2 gi + l 1 



^30 
°03 



4 



;,0 



(66) 



r(i)r ( ?*±l) + E! 



■l/2)r((2 gi -l)/4) 



r(3/2) 



(5 



(i)\2 



r(-(2g, + l)/4)r(-(2g t -3)/4) 
r((2 ft + 5)/4) 



(i)\2 



(i)\2 



One can check that the first term in the right hand side, is identical to l q ._i(0) as it should. 
This result is interesting in view of the crude approximations we had been relying on in 



[1,2]: In particular, the term of exp[— (b/x)\ of (60) was expanded at leading linear order in 
b ~ g(f{b) <C 1. One can see now that, were it not for the terms at = 0, this expansion 



yielded the correct leading order estimate of (44). 



One has now to deal with the remaining integrations, as they show up in (28), that is, 
firstly, with the subsequent integrations of 



+°° . a \ vf 

dV{ e E ~p 



dYt e~ im ^ G 



03 



\g(p(b)m 2 VjVp 2 
128N C E 2 p 2 



2ft + l 1 







(67) 



where, while having kept the same symbol a\, a\ as in (28), the original n-component vectors 



«i, a 2 have been promoted to iV-component ones, defining a, := (1, 0, 0, 1) ® ai, for i = 1,2. 
Rather than integrating a Meijer's function expansion which doesn't terminate, one can take 



advantage of formula 20.5.(8) in Ref. [16], so as to obtain for (67) the exact result of 



f g<pQ>) \ 



2 4 

m 



a\a\ 43 I 8N C V ®\a\ J E 2 p 2 



1 1 

2' 2' ' 

2qi + l 1 n 

4 ' 2' ' 



(68) 



Then, in view of (|28|), two extra integrations remain to be carried out, 

r+00 

00 



47r ^ 32 / i ( pcp(b) \ 2 m 4 
X ofaf 43 \ 8N C V oi\a\ J E 2 p 2 



1 1 

2 , 2 ' w > w 

2g»+l 1 n 

4 ' 2' U 



(69) 
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that can be most conveniently re-written as 



/+oo /"+OO 
da\e ia ^ 2 / daje^ 71 
oo i-oo 



/2 



341 i V^( & ) / ™ 4 



3-2 qi 1 i 

4 ' 2> ' 

i ^11 

2 ' 2 ' ' 



where use has been made of the inversion formula 5.3.1(9) of [15] . 



(70) 



I x- 1 



G%- I x 



1-6. 
1 - a r (j) 



(71) 



By symmetry, (|70|) reads indeed as 

30 

— Itt I da\ 



sin(ai7^/2) 



a 



~ da ,_sinKT72) 

— OO 



a 2 



(72) 



X 



23 /8iV, fa\aj\ 2 EY 
34 1 z \gip(b)J m 4 



3-2 9i 1 -i 

4 ' 2' ' 

i ^11 

2' 2' ' 



so that, by expanding the sme-functions and recalling the relevant properties of the T % - 
generators: Vi = 1, . . . , N, 3a = 1, . . . , n, such that (V) 2 = (T a ) 2 , (V) 2k = (T a ) 2k = (T a ) 2 , 
Vfc > 1, and Vo, Tr (T a ) 2 = C A = N c , one obtains 



-4**./ i]W in(a ' ,/2) 

oo \ j=1 



a 



(73) 



( a\ai \ 2 E 2 p 2 



3-2^ 1 1 
4 ' 2' ' 

1 1 



2 ' 2 ' ' 



i \g<p(b) J m 4 
where the trace over the internal color degrees has been taken. 



Recalling the shorthand Dn := N and getting back to (57), integrated over V{ 2 i one 
obtains therefore a 4-point fermionic amplitude of 



N 



64vr 



9\ - 
2 \ 4 



. JV , . JV* JV M 



(74) 



l<j<N 

x e ( ±x ) n [i 

monomials ^2lj = N(N—l) 



da\ 



sin or 



a 



j sina 2 C 2sl^[c ( ajaj \ 2 E 2 p 2 



da 3 2 



3 G 34 



z \g(f(b) J m 4 



3-2qj 1 1 

4 ' 2' ' 

i ^11 

2 ' 2 ' ' 
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In the second line, the sum runs over the monomials composing the Vandermonde determi- 
nant, and the (±1) is here to mean that each monomial appears with a plus or minus sign. 
Already at N = 4, one gets 729 monomials of overall degree 12 in the the eigenvalues £j 
that, after some cancelations, reduce to 201 terms. But at N = 6, even after reduction, one 
gets 56183 monomials of overall degree 30 that require several pages to be written down. 
Needless to insist on the number of monomials entailed in the Vandermonde determinant at 
N — 12, relevant to the case N c = 2, not to speak of N c = 3 for which N = 32. 



Before proceeding to the full integration of (74), it is worth observing that the last line, 

„2 



-sinx-G^/V^- 



(75) 



was the essential part of our first and crudely approximated estimate in Ref. pQ, and the 
result was 



4 

gy Jo 



fir 



dx 



sm ix) 



x 



(76) 



Based on (76), we had found for SU(2) that at large values of g<p(b), corresponding to small 
impact parameters, I(g<p) ~ 1/flV, obviously damped in comparison to the large impact 
parameter limit of gip <C 1, where I{g<p) ~ C st . The interpretation proposed in Ref. pQ 



and based on (76), preserves its essential features within the present paper's more accurate 
calculation where one gets a behavior of l/gtp at large impact parameter, and a behavior of 
^fgtp at small impact parameter. 

Focusing on a given integration, that of index j for example, one gets for the integrals of 



(74), 



sm 



dx 



x 



VSNc Ep 



X 



sm 



dy 



gtpVi m 2 
s/8N^ Ep 



Gil I x z y 



2„.2 



y 




(77) 



with {■ ■ ■ a r (j) ■ ■ ■ } and {• • • b s • ■ • } the sequences of numbers such as displayed in (74) 



These two sequences of numbers are such that the reduction formula 5.3.1.(7) of [15], can 



be used to simplify somewhat the Meijer's special function entering (77) and write 



3-2 qi 1 
4 '2' 
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2 ' 
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2 ■ 
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1, 1 
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3-2 9i \ 


) =Gg| 


x y 


2' 
1 

2' 
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(78) 
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This allows one to re- write (77) as 



dx sin 



gip\/i m 2 



dy sin 



VW>Ep 
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o, 




Gil 


x y 








o, 



(79) 



where the following identity, 5.3.1.(8) of [13], has been used, 



G pq | x 




(80) 



In the end, parameters are such, in the sine function as well as in the Meijer special 
function, the sequences of numbers {■ ■ ■ a r (j)+a • ■ ■ } and { 



b s + a ■■■}, that 20.5.(7) [IS], 



can be used twice so as to yield for (|79|) the result 

2 



gip(b)\/i m 2 
gcp{b)Vi m2 " 36 \ I v / 128iVc Ep 



7ry/8N c Ep 



1, 1, 



i i i 2^+3 1 1 
L i L i L i 4 ' 2' 2 



(81) 



where the inversion formula 5.3.1(9) of [TS], that is (71), has been used again in order to 
express the result in a form suited to a g(p(b) <C 1 limit expansion of the net amplitude. 
Gathering all factors, the overall omitted multiplicative constants included, the amplitude 



under consideration, (|74|) can eventually be displayed as 
J\f 



1 



£ 

monomials 



m 
s 



„ N N 

2 \ /Q4 n 2\ 4 



z)4 (-8i) — N r — +1 



(82) 



X< - T 36 



l<j<iV 

±1) n [l-ic-OT 

E ?j =jv(iv-i) 

^ \ 2. 

g(p(b)\/i m 2 /'s 



g<p(b)y/ij 



v/128AT c ^/l - m 2/s 



1, 1, 
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2 

2<ft+3 1 1 
)-"->-"-) 4 ' 2' 2 



1, 1, 1 



where a notation typical of perturbative Drell-Yan scattering processes has been introduced: 
's = (Pi + P2) 2 , the quark's 4-momenta, the Pi, being related to their respective hadronic 
carriers in the usual way, Pi = x{Pi, with Xi £ [0,1]. Beyond the notation itself, this is 
certainly far from an innocuous fact, as an explicit link shows up between a perturbative 
partonic content and a hadronic non-perturbative amplitude |18| . within one and the same 
overall expression. 

It is interesting to notice that the small g</>(5)-limit expansion of the G 42 ,-Meijer's functions 
is the more relevant as their argument is not gtpib) itself, but g(p(b) x (l/vT28AQ x m 2 /s, 



20 



that is a particularly small number in the eikonal limit. At leading order, corresponding 



to the value 1/2 of the 6 s -parameters, the two last terms of (82) combine into a constant, 
gy?(&)-independent contribution, as it should. Also, the combinations of initial parameters 
that have been generated througout the calculation are worth observing: They turn out to 
be the natural 'words' in term of which the amplitude is finally written. 

The concise form obtained for the fermionic correlation functions can certainly be ex- 
ploited in various situations, limits, and along various numerical and/or semi-analytical 
ways. Similarly, when the full flavor and mass and structure function richness of QCD is 
restored, these generic forms of amplitudes could be used and tested at a phenomenological 
level to describe hadronic collisions at intermediate energies, for example. This however, at 
the moment, falls beyond the scope of the present article where one is interested in under- 
standing quark's binding potential such as derived in [2]. As advertised in the main text 



of the present article, that result, basically given in (27), displays the occurrence of an in- 
triguing ^-parameter whose non-zero value is essential for obtaining a non-trivial confining 
potential. 



As the derivation of (27) relied upon frequently crude approximations, a natural issue is 
to know wether that ^-parameter is an avatar of the approximate calculational scheme used 
in [2], or not. It is this question that can be addressed now that, with respect to our first 
crude estimations [H [2] , the calculational accuracy has been substantially improved. 

Following the standard procedure that relates a 2-body eikonal amplitude to a potential 



acting between them [21 [13112], it turns out that from (82), the next-to-leading order of the 
g^p(b) <C 1 limit expansion is required. 



Now, as made obvious by Meijer's functions small g<f(b) expansions, ( ]63| ) and (65), the 
next-to-leading contributions are obtained from the parameters 2qj + 3/4 whenever qj = 0, 
and instead of order gtp(b), as derived previously [2], are indeed on the order of ^Jgip{b). 

However, as argued in [2], it is the logarithm of that next-to-leading order contribution 



which is associated with the calculation of V(r). In view of (17) though, one has (f(b) ~ 
e -\m ( see a i so JTT] , where a form of ip(b) ~ e~ b2 / Aa is proposed). Therefore, not solely that 
power of 1/2, instead of 1, but any other possible next-to-leading power of gtp(b) would ever 
bring the slightest change to the conclusion of [2]: A non-trivial quark confining potential 
requires a no n- vanishing ^-parameter value. 

In eikonal and quenching approximations at least, this conclusion is now reached on 
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a more accurate basis, and it is made accessible thanks to the non-perturbative Effective 
Locality property of QCD. 



V. DISCUSSION 

In eikonal and quenching approximations, the QCD property of effective locality has 
been exploited in this article to derive some interesting results. In the first place, evaluating 
QCD amplitudes, it has been argued that an awkward term of 5^ (6), appearing as an 
overall multiplicative factor in an exponential argument, could not be avoided and is in no 
way related to any approximation scheme. 

In a previous analysis [2], the position had been adopted that this term be regarded as 
some relic of an implicit perturbative input, somewhat unavoidable but irrelevant to the 
non-perturbative context [5] of the Effective Locality consequences. 

This term of 5^{b) was accordingly corrected and in a first attempt, turned into the 
smoother form of a gaussian distribution of the scattering quarks transverse distance, i.e., 
their collision's impact parameter. We could note that in an analysis of "Relativistic Quark 
Models in the Quasi-potential Approach" , a similar form had been proposed [17]. Against all 
odds, though, the quark confining potential associated to the scattering eikonal amplitude 
is then found to be identically zero [2]. 

An elementary deformation of the random transverse quark gaussian distribution has 
been suggested in order to prevent such an odd result from occurring [2]. Fortunately 
enough, that deformation (which, more precisely, doesn't bear on the Levy flight distribution 
itself, but on its characteristic function), falls into a canonical generalization of random 
variable gaussian distributions that are compatible with a generalized central limit theorem 
of statistical physics. 



Another interesting aspect of the amplitude's general form of Eq. (82) is that, within 
one and the same expression, it is able to encode the hadronic non-perturbative amplitude 
together with its perturbative partonic content. It is worth noting that such a link has been 
advocated recently (and the way by which this link appears) within the AdS/QCD approach 
to non-perturbative QCD [18J. 

Conclusions are here obtained on the bases of an eikonal and quenching approximation to 
the full QCD fermionic scattering amplitudes. This means that so long as quenching is not 
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relaxed, the relevance of that ^-parameter remains somewhat uncertain; and this, in spite of 
its obvious relation to color confinement compared to a possible transverse, non-commutative 
geometrical aspect of non-perturbative QCD. 

On the other hand, the results of [20] are not restricted to any sort of approximation. In 
this paper, using Random Matrices, it has been shown that the QCD fermionic amplitudes 
comply with that conjecture, and this might as well suggest that our present conclusions 
could be independent of the quenching approximation. Still, that remains to be checked. 

In these perspectives, quark propagator's full calculation and self energy effects should 
bring some new understanding, and will be the matter of forthcoming analyses. 



Appendix A: Proof of Eq. (62) 



Starting from I p (b) such as given in (60) and evaluated thanks to [14J, at p > and b > 0, 
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one proceeds to calculate (59), that is the integral 
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Now the first integral of (A4) can be evaluated by means of Cauchy's theorem: 



dz z v e 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



where the integration contour Y is as depicted in Fig. [TJ and where, to begin with, p must 
be restricted to positive rationales. In this way, zero is obtained as 



i P {c)-vr +1 i dxx>'< ■ 



0. 



(A6) 
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FIG. 1: Integration contour. 



Now, the analyticity of Meijer's special functions in their arguments, as well as the possibility 
of continuing their parameters to arbitrary complex values [15] allows one to define the 
required integral as 



_i _• 2 AW/yjg; 
dx x q% 2 e lx 



o 



and (A2) eventually under the form of 
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that is Eq. (62) in the main text. 
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